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Abstract 

We have developed a Green's function formalism to compute the local field 
distribution near an interface separating two media of different dielectric con- 
stants. The Maxwell's equations are converted into a surface integral equa- 
tion; thus it greatly simplifies the solutions and yields accurate results for 
interfaces of arbitrary shape. The integral equation is solved and the local 
field distribution is obtained for a periodic interface. 
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The investigation of the electric field distribution in heterogeneous dielectric media has 
been an established topic However, most studies have been restricted to interfaces of 
simple geometry. These studies are important because the detailed form of the local field 
distribution is responsible for the shifting of electronic states at semiconductor interfaces 
0. The same form is also important for interpretation of scanning electrostatic force mi- 
croscopy experiments However, real surfaces must have some degree of roughness. 
Such an imperfection not only influences the local field distribution but it also determines 
the implications on physical systems ||. 

Theoretical calculation of the local field distribution requires a detailed solution of 
Maxwell's equations. It is a formidable task for interfaces of complex geometry because 
of the complicated boundary conditions. Thus various attempts have been made in the 
literature to avoid the match of boundary conditions at each interface of the system ||. In 
this paper, we develop a Green's function formalism to compute the local field distribution 
for an arbitrary interface separating two media of different dielectric constants. Moreover, 
the local field distribution is calculated near a periodic interface. 

The object of this investigation is two-fold. Firstly, we aim at developing practical 
computational methods for calculating local field distribution near interfaces of arbitrary 
shape. Secondly, the impact of geometric resonance on local field distribution will be studied. 
The results show that the interface profile will provide an extra degree of freedom for tunable 
resonance enhancement of local field. 

Integral equation formalism. The electrostatic potential satisfies the Laplace's equation: 

V-[e(r)V0(r)] = -47rp(r), (1) 

with standard boundary conditions on the interface, where p(r) is the free charge density, 
e(r) equals €2 in the host and e\ in the embedding medium. 

In a conference paper [[/J, we developed an integral equation approach to compute the 
Green's function for an interface separating two media of different dielectric constants. By 
using the Green's identity, we were able to convert the volume integral equation into a surface 
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integral equation; thus it greatly simplifies numerical solutions and yields accurate results 
for interfaces of arbitrary shape. Here we re-iterate the formalism to establish notations. 

Let V\ and V2 be the volume of the embedding and host medium, separated by an 
interface S. Denoting 9(r) — 1 if r G V\ and otherwise, leads to an integral equation 0: 

[1 - ^(r)]0(r) = o (r) + £ jf dS' [n' ■ V'G(r, r')] 0(r'), (2) 

where u = 1 — e.\fe<2, n' is unit normal to S, G(r,r') = l/|r — r'| and </>o is the solution of 
V 2 0o = -47rp/e 2 . 

Accordingly, our approach aims to solve a surface integral equation for the potential at 
the expense of a two-step solution: 

1. step 1: determine </>(r) for all r G S by solving Eq.(^), and then 

2. step 2: obtain 0(r) for all r by using Eq.(g) and the results of step 1. 

In step 1, we encounter a singularity when the integration variable r' approaches the point 
of observation r. To circumvent the problem, we take an infinitesimal volume around r and 
perform the surface integral analytically, we find 

(l - |) 0(r) = o (r) + £- £ dS' [n' • V'G(r, r')] 0(r'), r G S, (3) 

where "prime" denotes a restricted integration which excludes r' = r. The (surface) integral 
equation (|3|) can be solved for 0(r G 5 1 ). 

Periodic corrugated interfaces. Here we apply the integral equation formalism to a peri- 
odic interface. Suppose the interface profile depends only on x, described by y — f(x), where 
f(x) is a periodic function of x with period L: f(x + L) = f(x). Thus medium 1 occupies 
the space y < f(x) while medium 2 occupies the space y > f(x) separated by the interface 
at y = f(x). The external field is E and 0o( r ) — — E • r. For a periodic system, 0(r) is 
a periodic function of the lattice vector T. In what follows, we adopt similar treatment as 
the Korringa, Kohn and Rostoker (KKR) method and rewrite the integral equation as: 

(l - |) 0(r) = -E ■ r + £- £ dS'G(r, r')0(r'), (4) 
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where the integration is performed within a unit cell. The structure Green's function (Gree- 
nian) is given by M: 



G(r, r') = V'G(r, r' + T) = £ g^L_i! ^1 . (5 ) 

^ ^ |r — r — mx| 2 



In view of the complicated calculus, the evaluation of the Greenian is a formidable task. 
However, by invoking the complex notation z = x + iy, z' = x' + if(x f ), we are able to 
evaluate the Greenian analytically: 

~ , , ^ 2idz'/dx' 27c[f'(x')sm27T(x-x')-smh27T(y-f(x'))] 

G(x,y;x , f (x )) = Re > = — 1 ; : ; ; —, — . 6 

v ' y ' ,jy " ^z-z'-m cos27r(x-a; / ) - cosh27r(y - /(x')) 

Eq.(^|) is readily derived by using the identity || for the infinite sum: 

OO J 

^ = 7TCOt(7Tz), 

z being a complex number, and by taking the real part of the result. Eq.(|6]) is a truly 
remarkable result - the analytic expression is valid for an arbitrary interface profile. For a 
planar interface, f(x) = f'{x) = 0, then G = and = on the interface, as expected. If 
the point of observation (x, y) is located at the interface, the Greenian has a finite limit as 



x' — > x: 



G(x, /(*); x', f(x')\x' - *) = 1+ / ( ^ ))2 . (7) 

We first solve Eq.(f|) for the potential 4>(x, f(x)) right at the interface: 

(l - |) <f>(x, f(x)) = -E f(x) + ±J dx'G(x, f(x); x', fix')). (8) 

Then we use Eq. @ to find the potential at any arbitrary point (x, y) S, using the potential 
at the interface. 

[1 - udifix) - y)]<p(x, y) = -E Q y + — J dx'G(x, y- x\ /(x , ))0(x' ) fix')), (9) 

where = for x < 0, and 9{x) = 1 for x > 0. The electric field is computed by the 
negative gradient of the potential. We find the electric field in medium 2: 



g „ = -g^l = Ee - £ /^ afi(j -"^- /W) ^,/W), (11) 



and similar expressions in medium 1. The partial derivatives of the Greenian can be com- 
puted analytically from Eq.(||). 

Solution by mode expansion. To solve the integral equation, we express the potential at 
an arbitrary point into a mode expansion: 

<f>(x,y) = EQ^Wa(2/ - /(*)), (12) 
jk 

where ipj(x) and are mode functions. The potential on the interface suffices: 

<Kx,f(x)) = EQ*W6(0) = 'EAtfjix), (13) 

jk j 

where Aj = J2k Cjk£k(fi)- Here a few remarks are in order regarding the choice of the mode 
functions. In theory, the choice of the mode function is somewhat arbitrary. In practice, 
these functions should be simple and easy to use. Common choice ranges from extended 
mode functions like the Fourier series expansions to localized mode functions like the step and 
triangular functions. Substituting the mode expansion Eq.flTB]) into Eq.(|j), the coefficients 
Ai satisfy the matrix equation: 



2J 4tt . 



A = -E V, (14) 



where 



Bij = J dxifji(x)ifjj(x), (15) 
M lj = J J dxdx'ij t (x)G(xJ(xy,x'J(x'))^(x') } (16) 
Vi = f dxiPi(x)f(x). (17) 



It should be remarked that the mode functions need not be orthonormal and the matrix B 
is non-diagonal in general. 
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Numerical results. As an illustration, we adopt the interface profile: 

2nx Anx 

fix) = — cos 0.1 sin , 

y 1 L L 

where the sine function is added to upset the reflection symmetry about x — 0. We adopt 

the step functions for the mode expansions: 

h h 

9i(x) = 1, Xi < x < Xi H — and 6^(x) = 0, otherwise, 

where h is the width of the step function. In what follows, we adopt 100 step functions, 
equally spaced in the unit interval x G (— L/2,L/2\. The integrals Eqs.(15)-(17) can be 
readily performed. 

We have attempted the following calculations: 

1. For positive ei/e 2 ratio, we find (f>(x, f(x)) right at the interface and compute S(f>(x) = 
4>(x, f(x)) — (f)o(x, f(x)). The planar interface result is 5(j) = 0. 

2. For negative ei/e 2 ratio, there exists a surface plasmon resonance. We find 8<j)(x) right 
at the interface near a geometric resonance of the interface. 

The geometric resonance is characterized by nontrivial solution of Eq.(14) even in the absence 
of Eq. Without loss of generality, we choose E = 1. Fig.l shows the results for 5<j>(x) and 
E y right at the interface. Two parameters of u — —1,-9 are used in the calculations. It is 
clear that the amplitude of S(f> increases with the dielectric contrast. It is more interesting to 
compute the potential near a geometric resonance of the interface. Fig. 2 shows the results 
for 8(f>(x) at a positive value of u — 3.3645. There is a giant enhancement of the potential. 
We also compute the y-component of the electric field at a fixed distance y = 1.1. For 
negative values of the ei/e 2 ratio, there is strong enhancement of the local electric field. In 
both cases, E y exhibits similar behavior as that of 5(f). It should be remarked that 5<f) and 
E y here are not computed at the same position. 

We have performed similar calculations for the triangular function and obtained essen- 
tially the same results. However, for the triangular functions, the B matrix is tridiagonal. 
Alternatively, we repeat the calculation by using the Fourier series expansion. 
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<i>(x,f(x)) = a + Y, 



2nnx 



+ b n sin 



2nnx 



) 



(18) 



a„ cos 



L 



L 



n 



We have adopted n — 1,2,3 harmonics in the Fourier series expansions. For u = —1, —9, 
the results (not shown here) of S(f>(x) agree with the calculations by using the step function 
expansions. However, the u = 3.3645 results exhibit some deviations. This can be clearly 
seen from Fig. 2 that the local field distribution contains rich harmonics near a geometric 
resonance of the interface. 

The Fourier coefficients for a few lowest harmonics are listed in Table I. From these 
coefficients, we reproduced 5(f) at x = and computed the error. As is evident from Table 
I, the error is small for the u — — 1, —9 cases, while there is a large error for the resonant 
case. Thus, to achieve more accurate results, higher harmonics need to be included in the 
resonant case. 

In summary, we have developed a Green's function formalism to calculate the local field 
distribution of an interface of arbitrary shape. The solution is accomplished first by solving 
the potential at the interface, followed by determining the potential at any arbitrary point 
by performing an integral over the interface. Thus it greatly simplifies and yields solutions 
for interfaces of arbitrary shape. Moreover, the local field distribution have been calculated 
near a periodic interface. 

This work was supported by the Research Grants Council of the Hong Kong SAR Gov- 
ernment under grant CUHK 4284/00P. K. W. Y. acknowledges initial collaboration with 
Professor Hong Sun. 
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Table I: Fourier coefficients for the potential at interface. 



constant cos27r:r cos47r:r cos67r:r sin27ra; sin47r:r sin67r:r error 

a ai a 2 a 3 b x b 2 b 3 

interface profile 0.0000 -1.0000 0.0000 0.0000 0.0000 -0.1000 0.0000 0.0000 

u = -1.0 -0.1739 -0.3268 -0.0337 0.0018 -0.0083 -0.0321 -0.0066 0.0004 

M = _9.0 -0.1228 -0.8065 -0.0246 0.0033 -0.0060 -0.0797 -0.0050 -0.0001 

u = 3.3645 12789 -6935 3822 -3005 2392 483 2676 -166 



9 



FIGURES 

FIG. 1. 5<f>(x) (thin lines) and E y (thick lines) plotted against x for two values of dielectric 
constrasts u = —1.0 and —9.0. 

FIG. 2. 6<f)(x) (solid line) and E y (dashed line) plotted against x for u = 3.3645. There is a 
giant enhancement of the local field. 
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